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A8STRACT 


The  theoretical  analysis  of  delta  wings  in  incompressible  flow  is 
treated  by  a numerical  1 i f * ing-surface  theory  based  upon  a velocity  po- 
tential formulation.  This  theory  has  been  expanded  to  include  spanwise 
velocity  effects  and  the  leading-edge  separation  associated  with  delta 
wings.  The  numerical  technique  has  made  use  of  both  fixed  and  free  wakes. 
In  the  fixed-wake  model,  the  analytical  results  of  Brown  and  Michael  have 
been  used  to  position  the  leading-edge  vortex.  In  the  more  refined  free- 
wake  model  the  leading-edge  separation  was  modeled  as  a discrete  number 
of  vortices  attached  to  the  leading-edge,  which  are  allowed  to  align  them- 
selves with  streamlines.  The  iteratively  determined  position  of  these 
vortices  resembles  the  experi mental ly-observ<  i spiral  form,  while  computed 
lift  coefficients  reflect  the  added  vortex  lift. 

Experimental  research  included  wind  tunnel  tests  to  determine  the  ver- 
tex location,  vortex  burst  location,  a" d upper  surface  pressure  distribu- 
tion on  a delta  wing,  AR  = 2,  in  an  oscillatory  airstream.  In  steady  flow 
an  axisymmetric  'tulip'  type  burst  was  observed  at  all  angles  of  attack; 
while  in  unsteady  flow,  the  vortex  core  gradually  dissipated  into  turbu- 
lence for  increasing  angles  of  attack  but  formed  a 'helical'  type  of  burst 
for  decreasing  angles.  The  pressure  distributions  showed  that  in  unstead, 
flow,  there  were  considerable  Dhase  lags  in  the  suction  peaks  associated 
with  the  vortices  but  virtually  no  lags  elsewhem  on  the  wing. 
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nomenclature: 


- i.'tin 

- area  of  vortex  box 

- aspe:t  ratio 
b - w i r j ^pan 

c - contour  variable 
: - rean  aerodynamic  chord  (-2/3  cp) 

-b  - root  ^hord 

C - aerodynamic  influence  coefficient 

- lift  :oe*ficient 

L 

2 x,  - normal  force  coefficient 

D - flow  field  domain 

k - reduced  frequency  (=-*—-) 

v - doublet  strength 

■ip  - potential  lift  constant 

<;  - vortex  lift  constant 

r - surface  normal 

- field  point 

Q - boundary  surface  point 
r - distance  between  P and  3 
s - local  winq  semi -span 
S - boundary  surface  variable 
t - time 
V - velocity 

WK AL  PH  - ratio  of  wake  inclination  angle  to  angle  or  attack 
( 3 i , S^,  S3)  - velocity  oriented,  normal  surfaim  axes 


IX 


( u,  v,  w)  - Cartesian  velocities 

(u,  uQ,  ur)  - polar  velocities 

(x,  y,  z)  - Cartesian  coordinates 

(x,  0,  P.)  - polar  coordinates 

(Xy,  Yy,  Zy)  - coordinates  of  vortex  core  endpoint 

iCp  - loading  coefficient 

Ap  - pressure  differential  across  surface 

a - angle  of  attack 

e - wing  semi -apex  angle 

> - vorticity 

9 - perturbation  velocity  potential 

* - angular  velocity 

Subscripts : 

i - lower  surface 
n - normal 
u - upper  surface 
s - surface 

(U,V,W)  influence  coefficient  subscripts  - influence  coefficients  of 
(u,v,w)  velocities 
W - wing 

V (matrix  subscripts)  - LE  vortex  shepc 
TE  - TE  vortex  sheet 
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INTRODUCTION 

Since  1974  the  Office  of  Naval  Research  (ONR)  has  been  conducting 
an  ir-deDth  study  of  vortex  flows.  The  aim  of  this  stud/  is  the  under- 
standing of  the  interaction  between  wings  and  vortex  'low  systems, 
f rom  the  understanding  will  come  advances  in  stall  and  buffet  bound? '’es. 
and  increased  performance  for  VTOL  and  STOl  craft.  Texas  'Jniver;  t.. 
is  under  contract  to  ONR  to  investigate  theoretically  and  experimento’ ' v 
the  ,C1  Ow  over  delta  wings.  This  report  includes  a complete  account 
: * the  researcn  bindings. 

- delta  wing  is  considered  the  starting  point  in  'r vesti aati ng  tne 
DOten*i.-.l  benefits  of  leading-edge  separation  because  of  its  simple  r’en- 
<fcm  and  because  the  vortex  flow  over  the  wirg  is  created  by  sel  f- induced 
seoarat;on.  In  addition,  abundant  experiment  evidence  is  available  to 
evaluate  theoretical  predictions  and  to  indicate  aspects  requiring  special 
attention,  theoretically  or  experimental ly. 

■"ne  available  experimental  evidence  leads  one  to  believe  the  flow  over 
a delta  wing  is  as  shown  in  Figure  1.  The  upper  and  lower  surface  boundary 
layer-s  separate  the  wing  along  each  leading  edge  to  form  fre^  shear 

layers.  These  shear  layers  then  roll  up  into  spiral  vortices  lying  above  ar.g 
inboard  c*  the  leading  edges.  The  suction  peak  on  the  wing  below  tn*>  pnmar. 
vortex  leads  to  seconcan  and  perhaps  tertiary  separations.  The  vorte*- 
inducec  suction  results  ii  increased  lift.  The  vortex-induced  lift  is 
nonlinear  with  respect  to  angle  of  attack  and  ^s  not  readily  predicted  b . 
linear  potential  flow  theories. 


( 

J 
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A complete  suwnary  of  the  theoretical  and  experimental  investigations 

14 

into  the  flow  about  delta  wings  has  been  compiled  by  Matoi. 

Correct  modelling  of  the  flow  should  include  the  effects  of  the 
primary  vortex,  secondary  vortex,  and  viscous  effects.  However,  expen- 
"iertal  evidence  [some  of  wnich  has  been  collected  at  Texas  AIH  University) 
ana  tne  success  of  theories  whicn  neglect  all  but  the  effect  of  the  primary 
vortex  suggest  cnly  the  first  phenomenon  need  be  included  in  the  theoretical 
model . 

Tne  theoretical  investigation  into  the  effect  of  leading-edge  separation 

? 

on  delta  wings  spans  over  a quarter  centu-'v.  In  1946,  2.  T.  Jones  , using 
slender  wing  theory,  postulated  that  low  aspect  ratio  wings  would  nave  a 
1 ' f t slope  of  - Aft/ 2 . This  value  is  only  slightly  conservative  at  low  angles 
of  attack  tut  since  the  theory  neglects  vortex  lift,  predicted  lift  falls 
increasingly  below  expenme'tal  lift  as  the  angle  of  attack  increases. 

Legendre^  in  1952  was  the  first  to  include  the  effect  of  leading-edge 
flow  separation  on  delta  wings  and  ascribe  the  significant  increase  in  lift- 
slope  to  this  phenomenon.  Using  conical  flow  theory  (as  do  most  analytic 
theories  developed  for  delta  wings)  Legendre  reduced  the  region  under 
cons  i He '•at  ion  to  the  cross-flow  plane  and  modeled  the  primary  vortex  sheets 


as  two  point  vortices  lying  above  the  wing.  The  position  and  strength 
of  these  vortices  were  determined  by  requiring  no  force  on  e’Ch  vortex 
and  a smooth  flow  off  the  leading  edge. 

4 

Brown  and  Michael  expanded  Legendre's  concept  by  adding  straight 
vortex  feeder  sheets  exte  ding  from,  the  leading  edge  to  the  isolated 


vortices,  as  seen  in  Figure  2.  These  sheets  accounted  .or  the  increaring 
strength  of  the  primary  vortices.  Brown  and  Michael  positioned  the  primary 


vortices  such  that  the  isolated  vortex  and  feeder  sheet  combined  were 


under  no  net  force. 


r /■  i 

Mangier  and  Smith3  and  later  Smith®  ar.u  Pullin'  allowed  each  feeder 
sheet  to  deform  into  a spiral  which  assym^totically  approaches  tne  vertex 
core*, Figure  3).  The  free-surface  condition  on  the  leading-edge  vortex 


sheets  was  more  rigorously  applied  by  separately  requiring  no  force  to  oe 
sustained  by  the  core  and  no  pressure  discontinuity  at  a finite  numDer  of 
points  on  the  sheet.  Smith's  work  is  especially  note-worthy  since  he 
cculd  easily  vary  tne  extent  of  the  leading-edge  sheet  outside  the  vortex 
core. 

D 

Nangia  and  Hancock  in  1968  were  the  first  to  apply  three-dimensional 
njrterical  1 i ft i ng-surfact  theory  to  delta  wings  using  the  model  of  Brown 
and  Michael  for  the  leading-edge  separation.  Their  procedure  involved  a 
collocation  method  for  determining  the  Four  fr  coefficients  for  the  assumed 
series  distribution  of  vorticity  over  the  wing.  Their  theory  did  include 
the  Kutta  condition  at  the  trailing  edge  and,  hence,  is  a departure  iron 
the  conical  flow  assumption  of  analytic  theories. 

Mook  and  Maddox^,  Kandil,  Mook  and  Nayfeh^  and  Rehbach^  represented 
the  leading-edge  separation  as  a system  of  discrete  vortices  wnich  are 
aligned  with  streaml  ines ( Figure  4).  Although  accurate  in  their  prediction 
of  lift  coefficients,  the  success  of  Kandil  et.  al.  at  reproducing  the 
experimentally  observed  vortex  snape  does  not  lend  confidence  to  the-r 
empirical  procedure. 

12 

Most  recently,  Weber,  et.  al . applied  a potential  flow  computational 
technique  to  the  problem,  -epresenting  the  wing,  rolled-up  vortex  sheets 
and  the  wake  by  piecewise  continuous,  quadraticclly-varying  strength  doublet 
panel  s ( "igure  5).  Their  technique  allows  the  tangent  flow  and  free-:.  jr;  - 
conditions  to  be  applied  independently  to  the  appropriate  panels,  and  hen.-. 


* tt-1 


satisfy  the  Kutta  condition  at  the  leading  and  trailing  edges. 


Finally,  Polhawus'^  concept  should  be  mentioned,  not  for  its  rigorous 
mathematical  justification,  but  because  of  the  accuracy  and  simplicity 
of  its  prediction  of  the  lift  of  delta  winqs.  The  concept  is  based  on  an 
analogy  between  the  vortex  lift  and  the  leading-edge  suction  associated  with 
the  potential  fiow  about  the  leading  edge.  Polhamus  defines  the  lift  of  a 
wing  with  leading-edge  separation  as: 

CL  = Kp  sina  os  ^ + Ky  sin^  COSa 

The  values  of  Kp  and  Ky  are  only  functions  of  aspect  ratio  and  can  be 
determined  by  linear  potential  flow  theory.  Graphs  of  them  for  deltt  wings 
are  giver,  in  Figure  6. 

Aerodynamic  data  on  wings  of  low  aspect  ratio  (of  order  1)  delta  wings 
with  sharp  leading  edges  is  well  documented. ^ 27 Date  on  wings  of  Aspect 
ratio  1.5  and  above  is  considerably  more  scarce  and  data  for  unsteady  flow 
past  these  wings  virtually  non  existent. 

The  steady  flow  is  characterized  by  two  vortices,  caused  by  leading 

edge  separation,  which  lie  above  the  wing  close  to  the  leading  edges.  With 

increasing  angle  of  attack  these  vortices  increase  in  strength  and  move 

inboard  causing  two  large  suction  peaks  in  the  spanmise  loading.  At  some 

point  along  the  length  of  the  vortices  the  phenomena  of  vortex  breakdown  i.e. 

the  brusting  of  the  tightly  rolled  vortex  core,  occurs.  At  moderate  angles 

of  attack  this  is  well  downstream  of  the  trailing  edge  of  the  wing  but 

with  increasing  angle  the  burst  point  moves  rapidly  forward  until  it  occurs 

(32-371 

above  the  wing  producing  unsteadiness  in  the  loading,'  ' Increasing  the 
aspect  ratio  causes  the  burst  to  reach  the  trailing  edge  at  low*r  angles  of 
attack. 
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Two  distinct  types  of  breakdown  have  been  observed1,3  ’ , one  more 

common  at  high  Reynolds  numbers  is  an  axisymetric  burst  where  the  core 
suddenly  expands  to  form  a "tulip"  shaped  "bubble".  In  the  other,  more 
generally  seen  at  low  Reynolds  numbers,  the  core  has  a sharp  kink  and  starts 
swirling  in  a helical  path  about  its  axis  before  breaking  down  into  turbulence. 

In  unsteady  flow  the  vortex  cores  move  both  laterally  and  vertically 
and  their  loci  have  been  measured  during  various  motions  for  low  aspect  ratio 
wings^88"^.  However  as  yet  no  effort  has  been  made  to  ascertain  the  effect 
of  flow  unsteadiness  on  vortex  bursting  or  to  determine  the  effect  of  the  burst 
on  the  transient  loading.  The  current  tests  were  planned  to  investigate  these 
areas . 
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I.  THEORETICAL  INVESTIGATION 
Aerodynamic  Theory 

Steady,  inviscid,  irrotational , and  incompressible  fluid  flow  Is 
characterized  by  a perturbation  velocity  potential  #(P)  satisfying  La- 
place's equation,  i.e.: 

v2$  = 0 (1) 

Green's  theorem  relates  the  value  of  ${P)  inside  a domain  D (I.e.  PcD) 
to  the  value  of  $ and  its  normal  derivative,  on  the  fluid  boundary  S: 

»0>>  = ,/  ^ <-&  dS  * ,/  «(«)  ± (^)  dS  (2) 

where  r is  the  scalar  distance  between  the  surface  point  Q and  P.  If  the 
boundary  S is  a two-dimensional  surface,  then  we  may  combine  the  integra- 
tion over  the  upper  and  lower  surfaces  as: 

^ $ e 1 

♦<p>  = <T5»* 


JL  (_L ) 

u 3nu  '4irr' 


+ $ -L_  (_L)  ds 

an.  Mwr' 

i 


(3) 


where  the  bar  over  the  surface  variable  indicates  integration  over  one 
side  of  the  boundary. 

But  n^  = -n^,  and  for  tangential  flow  to  occur  on  the  upper  and  lower 
surfaces : 


3n 

u 


(4) 


Then  the  velocity  potential  can  be  written  in  terms  of  its  upward  normal 


derivative  as: 


(5) 


s<p> - _'<♦ , • »«>  ir  O 

s u 

The  velocity  induced  at  a point  is  the  gradient  of  the  perturbation  velo- 
city potential,  v$(P).  Since  this  operation  is  performed  in  the  domain  D, 
the  operator  can  be  brought  inside  the  boundary  integration.  Introducing 
K as  the  discontinuity  in  potential  between  the  lower  and  upper  surfaces. 


K = i - $ 
t u 


then  the  induced  velocity  can  be  written: 


(6) 


*(P)  - - (f) 


(7) 


If  S is  divided  into  a number  of  small  subregions  ? , such  that  K can 


be  considered  constant  within  these  subregions  then: 

K. 


v<p)  ■ <-  / ' i * 


r) 


(8) 


The  surface  integral  over  each  subregion  has  been  shown  (Ref.  15)  to 
be  equivalent  to  a contour  integration  around  the  subregion.  Further, 
this  contour  integration  represents  the  influence  of  a vortex  filament  of 
unit  strength  on  the  contour  of  the  subregion,  as  illistrated  in  Figure  7. 
The  contour  integral, 


C.j  « - „ / c x v (7)  dc  = - / v ^7  (l)  dS 


c ■ 


■ n r 


(«)) 


J 5j 

becomes  the  vector  aerodyramic  influence  coefficient  of  subregion  j on  the 
point  P- . Using  the  subscript  i of  V to  denote  the  velocity  induced  at 
the  point  , the  total  induced  velocity  is  written: 

V.  - F K.C.  . 

1 4”  j J 


(10) 
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In  matrix  notation,  the  velocities  at  several  points  are: 

<Vt  = j;  [C]  4Kt 


(ID 


For  the  induced  velocity  field  given  by  (11)  to  be  accurate,  the  aerody- 
namic influence  matrix  must  adequately  represent  the  lifting  surface  by 
including  all  of  the  surface  over  whic*-  « discontinuity  in  potential  exists, 
sufficiently  approximate  the  geometry  of  the  surface,  and  be  consistent  with 
the  assumption  that  K is  constant  within  each  subregion. 

The  distribution  of  K over  the  surface  is  determined  by  three  con- 
ditions - undisturbed  flow  at  infinity,  tangent  flow  over  the  boundary, 
and  no  pressure  di'continuity  across  a free  surface.  The  condition  of  un- 
disturbed flow  at  infinity  is  implicit  in  the  Green's  theorem  derivation 
of  the  potential  equation.  It  can  also  be  interpreted  as  meaning  K is 
zero  upstream  of  the  leading  stagnation  point.  When  K is  so  specified, 
the  influence  coefficients  of  the  stream  surface  upstream  of  the  leading 
stagnation  point  need  not  ue  calculated. 

Over  the  wing  and  wake,  the  flow  must  be  tangent  to  the  lifting  sur- 
face, i.e.,  the  net  velocity  normal  to  the  surface  must  be  zero.  The  free- 
stream  velocity  component  normal  to  the  surface  at  point  Q..  is: 


(V  ).  = n.  • V 
v n«  i i 


and  the  induced  normal  velocity  at  point  Q.  is: 


(12) 


<Vi  ■ h L"i 


CjJ  (K)  ■ 


4* 


<K> 


(13) 


For  the  net  normal  velocity  to  be  zero,  the  induced  velocity  must  be  the 
negative  of  the  freestream  component,  i.e.. 


1C  J.  (K)  = -4it  (n.  • V ) 
L nJi  i °° 


(14) 


For  several  collocation  points  at  which  the  tangent  flow  condition  is  sat- 
isfied. 
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[C  ] { K}  *■  -4*  in  • ¥ k (15a) 

The  aerodynamic  influence  matrix,  [Cn3.  can  be  calculated  from  the 
known  geometry  of  the  lifting  surface,  while  the  frcestream  velocity  and 
surface  normals  are  also  specified.  Thus,  using  as  many  collocation  punts 
as  unknown  , the  doublet  distribution  over  the  surface  can  be  computed. 

Usually,  each  subregion  of  the  surface  is  provided  with  its  own  collo- 
cation point  located  near  the  centroid  of  the  subregion.  Empirical 
schemes, with  anc*  without  justification  have  been  used  to  strateqi 
cally  place  the  control  points  at  positions  other  than  the  centroid, 
lhese  alternate  positions  can  improve  a program’s  results,  since  the  lo- 
cation of  the  collocation  point  within  the  vor  .ox  box  can  effect  the  doub- 
let strength  of  the  box. 

The  third  condition  which  the  free  portions  of  the  lifting  surface 
must  meet  is  that  of  no  pressure  differential  across  the  surface.  The 
vorticity  on  the  surface  is  given  by  the  gradient  of  the  discontinui ty 
in  potential: 


y = n x (V£  - Vy)  = n x vK 
The  loauing  on  the  surface  is: 

aP  = pV  x y 


from  which  the  loading  coefficient  is 


= 2 


(n 


v)  *- 

' V 


(16) 


(17) 


(18a) 


(18b) 
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r . 2 , 3K  . u 3K  . M 3Kv 

nCp  - 77  <u  37  * v v " 57> 


(18c) 


If  a set  of  axes  (S^  S2.  S^)  is  placed  on  a free  doublet  surface  and 
oriented  with  the  velocity  vector  as  shown  in  Figure  8.  Then: 

V « |v|  s. 


n = S3 

'■37  *1  ♦ **  ♦ H7  S 


n x v = - S ♦ rf-  S. 

^ tv)  |^S3 


(188) 


For  a free  surface  AC_  = 0,  therefore, 

%*0  (19) 

From  (19),  it  is  apparent  that  on  a free  surface,  K is  constant  along  a 
streamline.  Thus,  in  steady  flow,  the  doublet  strength  at  a point  Q in 
the  wake  equals  the  doublet  strength  at  the  trailing  edge  point  upstream 
of  Q.  This  is  merely  a restatement  of  the  requirement  that  free  vorticity 
be  aligned  along  a streamline  in  steady  flow.  More  importantly,  it  also 
implies  that  if  two  sides  of  the  subregions  on  a free  surface  coincide  with 
streamlines,  the  subregions  degenerate  into  semi-infinite  strips.  This  is 
illustrated  in  Figure  9. 

Instead  of  providing  collocation  points  in  the  wake  and  using  the 
tangent  flow  condition  to  solve  for  the  K distribution  in  the  wake,  the 
influence  of  the  wake  strips  may  be  added  to  the  appropriate  edge  subre- 
gions, boxes,  (as  seen  in  Figure  10)  since  their  strengths  are  eoual.  This 


Si 


is  desirable  numerically  since  it  reduces  the  nuirber  of  simultaneous  equa- 
tions to  be  solved. 

Another  numerical  simplification  results  from  symmetry  conditions. 
Whenever  one  box  (the  image)  is  known  to  have  the  same  strength  as  another 
(primary)  box,  the  influence  of  the  image  box  can  be  added  to  the  influence 
of  the  primary  box.  This  simplification  can  result  from  symmetry  such  as 
that  occurring  between  right  and  left  wings  of  an  aircraft  in  straight, 
level  flight.  Such  symmetry  includes  the  tangent  flow  condition  whicn 
need  only  be  applied  to  the  primary  box. 

For  a high-aspect  ratio,  planar  wing  in  rectilinear  flight,  it  has 

been  found  necessary  to  rigidly  enforce  both  the  tangent-flow  and  free- 

1 8 

surface  conditions  on  the  t. ailing  wake.  In  such  cases,  wake  roll-up 
can  be  neglected  and  a planar  wake  consisting  f semi -infini te  strips 
trailing  back  in  the  x direction  can  be  assumed. 

Numerical  Applications  and  Results 
Attached  Flow  Model 

The  lifting-surface  theory  is  applied  numerically  to  delta  winqs  using 
the  sign  convention  shown  in  Figure  11.  The  wing  lies  in  the  x-y  plane 
with  the  x axis  pointing  aft  along  the  line  of  symmetry.  For  a wing  as 
described,  the  normal  vector  is  k which  makes  the  particular  equations  for 
tangent  flow  and  loading  coefficient: 

[Cj  • K 1 = -4”  iw^ - (15b) 
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If 

i The  wing  is  divided  into  subregions  (vortex  boxes)  by  chordwise  and 

} spanwise  lines  which  divide  the  interior  of  the  wing  into  rectangles  and 

form  triangles  at  the  leading  edge,  as  shown  in  Figure  12.  Collocation 

f 

points  lie  along  spanwise  lines,  halfway  between  the  left  and  right  sides 
of  the  boxes.  These  positions  which  are  the  centroids  of  the  rectangles, 
but  not  of  the  triangles  were  employed  to  facilitate  the  calculation  of 
the  gradient  of  the  doublet  strength.  Comparison  of  results  using  this 
scheme  and  one  using  the  centroids  of  the  leading-edge  triangles  as  col- 
location points  yielded  negligible  changes. 

By  adding  a wake  consisting  of  semi-infinite  extensions  of  the  chord- 
wise  lines,  an  attached  flow  (without  leading-edge  separation)  model  re- 
sulted. This  model  was  made  more  realistic  by  allowing  the  planar  trailing- 
edge  wake  to  incline  upwards  at  an  angle  described  in  terms  of  the  ratio 
of  wake  inclination  angle  to  angle  attack  (WKALPH) . For  values  of  this 
ratio  between  zero  and  unity,  negligible  changes  in  lift  resulted.  For 

the  results  to  be  presented  in  this  study,  this  ratio  has  been  set  at  one- 
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half  following  Gersten  and  Bollay  . 

The  computed  lift  slope  of  this  model  was  within  5%  of  that  calculated 
for  Polhamus'  potential  flow  model  when  an  8 x 8 grid  (36  boxes  on  each 
k half-wingj  was  used,  and  within  1%  when  a 10  x 10  grid  (55  boxes)  was  used. 

Thus,  a model  with  a 10  x 10  grid  has  been  used  for  the  reported  computa- 
tions. 

Two  models  of  the  leading-edge  separation  are  presented  here.  One  is 
a fixed-wake  model,  where  a simplified  form  of  the  primary  vortex  is  assumed 
and  the  pressure  distribution  is  calculated  without  checking  the  correct- 
ness of  the  assumption.  This  model  is  designated  the  Brown  and  Michael 
model  since  it  closely  matches  their  assumed  form  and  uses,  indirectly. 


their  analytical  calculations.  The  second  model  is,  to  a large  extent,  a 
free-wake  model.  That  is,  a guess  is  made  of  the  shape  and  strength  of 
the  leading-edge  vortex  sheet;  and  this  guess  is  iterated  to  improve  the 
leading-edge  sheet  until  it  satisfies  the  tangent-flow  and  free-surface 
conditions.  This  model  in  some  respects  matches  that  of  Mangier  and  Smith 
and  also  uses,  indirectly,  the  analytic  calculations  of  Smith. 

Calculation  of  Pressure  Distribution 

For  a planar  lifting  surface,  such  as  the  attached  flow  model,  the 
loading  at  any  point  is  dependent  on  3K/ax  and  the  lift  per  unit  span  is 
proportional  to  the  doublet  strength  at  the  trailing  edge.  This  is  true 
because  the  tangential  velocity  over  the  lifting  surface  is  constant.  The 
primary  vortex  sheets  of  a delta  wing,  though,  induce  significant  veloci- 
ties which  vary  over  the  wing.  For  each  point  at  which  the  loading  is  to 
be  computed,  the  two-dimensional  gradient  of  the  doublet  strength  and  the 
tangential  velocities  must  be  found  as  suggested  by  Equation  18e. 

At  present  the  gradient  is  calculated  using  three-point  finite  dif- 
ferences. These  differences  use  the  position  of  the  collocation  points 
and  their  associated  K values.  For  points  near  the  edges  of  the  wing, 
either  forward  or  backward  differences  are  used,  and  the  symmetry  condi- 
tion 3K/ay  = 0 is  incorporated  at  points  near  the  centerline. 

The  induced  velocities  are  calculated  using  the  influence  of  the 
leading-edge  sheet  only,  since  the  planar  wing  cannot  induce  tangential 
velocities  on  itself  and  the  influence  of  the  trai 1 ing-edge  wake  is  as- 
sumed to  be  small . 

The  normal  force  coefficient  is  the  integrated  loading  distribution, 
which  is  approximated  by: 


I ACpi  A< 

CN  = h-A  — 

i 1 


whe»v  A.  is  che  a.-®a  of  each  box. 


'he  lift  is  then: 


— CN  COSa 


Similarly,  the  moment  about  the  y axis  't: 

-E  iCpi  A.  x. 

r = 1 (22 

UMLE  I A.  Z v 

i 1 

where  the  location  of  the  collocation  point  is  used  for  the  moment  arm 
( x • ) of  each  box. 


Fixed-Wake  Model 

The  model  of  Brown  and  Michael  was  duplicated  by  placing  a series  of 
planar  boxes  along  the  leading  edge  which  extend  to  the  position  of  the 
vortex  core  predicted  analytically  by  Brown  and  Michael.  (See  Figure  13.) 
The  box  sides  extending  from  the  lead;ng  edge  are  perpendicular  to  it  and 
reprint  the  feeder  vortices  of  Brown  and  Michael.  The  sides  of  the  boxes 
opposite  the  leading  edge  form  two  concentrated  vortices  representing  the 
primary  vortex  system.  Experimental  observation  has  shown  that  the  core 
can  be  approximated  by  a straight  line  which  extends  from  the  apex  to  the 
trailing  edge  and  then  turns  toward  the  freestream  direction. 

In  accordance  with  the  treatment  of  the  Kutta  condition,  these  boxes 
have  the  same  doublet  strength  as  the  leading-edge  box  they  abut.  In  a 
strict  sense,  since  the  vortices  do  not  leave  the  edge  in  the  streamwise 
direction,  the  Kutta  condition  was  enforced  at  neither  the  leading  edge 
nor  the  trailing  edge,  where  the  attached- flow  model’s  wake  was  retained. 
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However,  the  numerical  model  had  enough  latitude  with  respect  to  the 
application  of  the  Kutta  condition  that  even  gross  approximations  concerning 
the  flow  as  it  leaves  the  leading  and  trailing  edges  yielded  acceptable 
pressune  distribution  near  the  edge. 

Tne  computational  procedure  is: 


1 . 

Input  AR,  _,  Xv,  Yy,  Zy 

2. 

Calculate  [CwJ 

3. 

Calculate  iw  : 

or 

4. 

Solve  Equation  15u  for  •. re- 

5. 

calculate  [C  ] and  [C ,] 
u v 

6. 

Find:  {u:  = 1/4-  [C  ] *K-  ♦ :u  t . j:  = l/4r 

[cv] 

7. 

Calculate  the  pressure  distribution  u ver  tne 

wing  using 

Equation  loe 

Resul 

ts  obtained  with  this  model  indicated  its  lirr.i 

tations  are 

the 

same  as  those  of  the  analytic  theory  of  Brown  and  Michael,  making  it  applies: le 
only  to  low  aspect  ratio  wings.  For  an  aspect  ratio  of  one.  the  predicted 
versus  curve  is  very  close  to  experimental  curves  (see  figure  14).  'hr 
agreement  drops  off  rapidly  with  increasing  aspect  ratio,  indicating  a pract.ca' 
limit  for  the  model  of  AR  * 1.5. 

Because  the  trai 1 i ng-edge  Kutta  condition  was  not  strictly  enforced,  the 
center  of  pressure  was  not  as  fa1"  forward  as  the  experimental  position.  As 
Figure  15  shows,  the  center  of  pressure  was  just  ^orward  of  the  centroid  o4  tne 
wing  whicn  is  the  center  of  pressure  predicted  by  slender  wing  theory. 
Interestingly,  wings  witn  practical  thicknesses  have  centers  of  pressure  closer 
to  those  predicted  by  the  1 i fti ng-surf ace  tneory  than  very  thin  wings. 


! 

l 


16 


The  computed  pressure  distribution  bears  a strong  resemblance  to 
experimental  results  as  illustrated  in  Figure  16.  This  figure,  showing 
tne  spanwise  pressure  distribution  near  the  center  of  pressure  (65X  root 
chord),  indicates  good  agreement  between  the  computed  and  measured  values 
of  the  centerline  loading  coefficient  UCp)  and  the  peak  aC,.  Because  the 
analytic  results  of  Brown  and  Michael  predict  a vortex  core  on  location 
outboard  of  experimental  positions,  *ne  suction  peak  computed  here  also  lies 
outboard  of  the  observed  peak.  Secondary  separation  will  also  contribute  to 
the  slight  discrepancy  between  the  two  curves  since  it  puts  a depression  on 
the  outboard  side  of  the  suction  peak. 

The  value  of  tne  fixed-wake  model  lies  in  its  simplicity.  The  required 
input  can  be  fitted  on  one  card,  ard  no  time-consuming  iteration  is  needed. 

In  fact,  the  run  times  were  less  then  5 seconds  (including  WATFIV  compilation) 
on  the  Amdahl  470V/6. 

Some  thought  has  been  given  to  improving  the  fixed-wake  model.  It  is 
believed  that  the  use  of  more  realistic  vortex  core  locations  would  move  the 
peak  suction  inboard  with  sub-equent  closer  agreement  with  experiment.  The 
aspect-ratio  range  for  which  the  program  is  accurate  might  be  extended  by 
adjusting  the  angle  at  which  the  feeder  vortices  leave  the  leading  edge, 
because  the  lift  coefficient  in  the  free-wake  model  was  found  to  be  very 
sensitive  to  this  angle.  This  adjustment  could  use  another  angle  for  all 
aspect  ratios  or  an  angle  which  is  a function  of  aspect  ratio. 

Free-Wake  Model 

Essentially,  the  representation  of  the  leading-edge  vortex  sheet  in  the 
■‘ree-wake  model  is  oy  means  of  discrete  vortices  attached  to  the  leading 
edge.  These  vortices  are  allowed  to  align  themselves  with  streamlines  to 
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meet  the  wake  criteria  of  tangent  flow  and  of  being  a free  surface.  These 
discrete  vortices  are  formed  o y the  boundaries  or  the  semi -infinite  subregions 
into  which  the  vortex  sheet  is  divided.  The  vortices  are  descrioed  by  line 
segments  whose  endpoints  are  specified  in  vertical  planes  which  coincide  w:t». 
the  spanwise  lines  that  bound  the  vortex  boxes  on  the  wing  surface  (Figure  17;. 
Beyond  the  trailing  edge  of  the  wing,  the  vortices  are  straight  lines  trail  115 
straigt  back. 

In  this  model,  the  trailing  wake  was  also  free  to  a certain  extent  to 
align  itself  with  streamlines.  Instead  of  being  a straignt  line  trailing 
aft,  as  in  the  attached-flow  model,  the  vortices  in  this  trai 1 ing-edge  wake 
consisted  of  a short  line  segment  attached  to  the  trailing  edge  and  tnen  a 
semi-infinite  trailer.  These  short  segments  extended  a distance  of  1C  of  tnt 
root  chord  aft  of  the  trailing  edge,  and  togw‘  'er  with  the  trailer  were  inclined 
vertical!.*  as  in  the  attached-flow  model  (Figure  18).  Within  this  inclined 
Diane,  trough,  the  short  segments  were  aligned  in  the  direction  of  the 
streamlines  leaving  the  trailing  edge.  Numerically,  no  problems  resulted 
fron  this  representation,  although  in  some  cases  the  short  segments  intersected, 
which  is  phy«ically  impossible. 

In  the  computer  program  the  doublet  strength  (K)  distribution  was 
not  obtained  dir?c:ly  through  Equation  15b.  Rather,  the  aerodynamic 
influence  matrix  was  broken  up  into  three  matrices  representing  the 
contributions  of  the  wing,  leading-edge  sheet,  and  trailing-edge  wake. 

Eauation  15b  was  rewritten: 

CUw  tK}W  = ' IVv  fK;V  ■ ^Vtl  K;TE  ;23) 

The  matrix  subscripts  refer  to  the  wing  (W),  leading-edge  vortex  sheet  (V), 
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and  tralling-edge  wake  (TE).  In  this  form  the  wing-on-wing  influence 
matrix,  [C^Jy,  which  does  not  change  with  each  iteration  and  which  is  the 
largest  matrix,  need  only  be  calculated  once.  This  matrix  can  then  be 
decomposed  by  the  Cholesky  matrix  solution  method  allowing  repeated,  rapid 
calculation  of  {K}y  for  each  different  right-hand  side  that  occurs  every  time 
the  wake  geometry  is  changed.  The  right-hand  side  calculations  are  performed 
using  the  assumed  geometries  and  strengths  for  the  leading-edge  and  trailing- 
edge  wakes. 

To  improve  the  assumed  geometry  of  the  wakes  the  velocities  at  the 
endpoints  of  the  segments  forming  the  wake  vortices  are  calculated  from 
the  free  stream  and  the  induced  velocities,  i.e., 

4*  {V}  = 4tt  {VJ  ♦ [C]w  {K)w  + [C]v  { K } y + [C]TE  {K)te  (24) 

These  velocities  are  then  used  to  improve  the  shape  of  the  wake,  so  that  the 
vortices  more  closely  follow  streamlines.  It  should  be  noted  that  the 
aerodynamic  influence  matrices  in  Equation  24  refer  to  the  influence  of  the 
lifting-surface  on  the  points  1r.  the  leading-edge  sheet.  The  matrices  in 
Equation  23  refer  to  the  influence  of  the  lifting-surface  on  the  wing. 

These  matrices  are  not  the  same. 

Originally,  the  updating  of  the  vortex  positions  in  the  leading-edge 
sheet  was  dependent  on  the  Cartesian  velocities  (u,v,w).  However,  for 
points  near  the  apex  or  close  to  the  vortex  core  the  curvature  of  the 
spiral  sheet  becomes  so  severe  that  the  Cartesian  velocities  are  useless 
in  predicting  the  streamlines  for  any  practical  grid  fineness. 

A great  deal  of  this  difficulty  is  overcome  by  adopting  the  polar 
coordinate  system,  shown  in  Figure  19,  which  is  consistent  with  the 
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experimentally  known  shape  of  the  vortex  sheet.  The  v and  w Cartesian 
velocities  are  transformed  to  tangential  and  radial  velocities  with 
respect  to  the  position  of  the  vortex  core.  The  variation  of  the  polar 
velocities  along  a vortex  is  much  smoother  than  the  variation  of  tne 
^artesian  velocities,  which  indicates  that  this  coordinate  system  is  more 
natural  for  predicting  the  vortex  geometry. 

Referring  to  Figure  £0,  the  new  angular  and  radial  coordinates  of 
the  downstream  end  of  a segment  are  predicted  with  these  equations: 


In  this  way  the  change  in  the  angular  and  radial  coordinates  between  any 
two  chord  stations  is  approximated  by  the  average  rate  of  change  of  each 
coordinate  multiplied  by  the  chordwise  distance. 

Unfortunately,  this  coordinate  system  does  not  completely  solve  the 
problem  of  increasing  curvature  near  the  vortex  core.  This  problem  is 
compounded  by  the  poor  prediction  of  the  velocities  near  the  vortex  core 
which  is  inherent  in  discrete  vortex  representations.  Consequently,  the 
sheet  near  the  center  is  poorly  reproduced  and  the  positions  of  the  vortices 
sometimes  indicate  an  unrealistic  intersection  of  the  sheet  upon  itself. 

To  circumvent  these  problems  a finite-sized  ortex  core  was  defined 
as  in  Figure  21.  In  the  program  two  circles  centered  about  the  vortex 
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core  limit  the  extent  of  the  vortex  sheet.  The  size  of  these  circles  grow 
continuously  downstream  of  the  apex  so  that  the  ratios  of  circle  radii  to 
local  wing  semi -span  is  constant.  If  the  endpoint  of  a segment  is  predicted 
to  lie  inside  the  inner  circle,  the  segment's  endpoint  is  moved  to  the 
center.  This  capture  is  reversible.  If  the  upstream  end  of  a segment 
(whose  downstream  end  is  captured)  lies  outside  the  outer  circle,  the 
position  of  its  downstream  end  will  be  approximated  such  that  the  downstream 
end  can  escape  the  core.  Additionally,  to  prevent  the  sheet  from  inter- 
secting itself,  a maximum  allowable  angle  is  specified,  as  also  illustrated 
in  Figure  21.  If  a vortex  exceeds  this  angle  it  will  be  moved  to  the 
center. 

5 

This  core  model  resembles  that  of  Mangier  and  Smith  , who  showed  that 
the  effect  of  the  core  of  the  spiral  sheet  on  the  field  outside  was  equi- 
valent to  a concentrated  vortex  at  tfe  center  of  the  core.  Smith6  later 
refined  their  work  including  the  specification  of  a maximum  sheet  extent 
angle  and  evaluated  the  effect  of  this  parameter. 

It  is  interesting  to  note  that  if  the  core  is  given  a large  enough 
radius  the  result  will  be  a model  similar  to  the  Brown  and  Michael  model. 

No  criteria  have  been  developed  here  to  predict  the  position  of  the 
vortex  core  due  to  the  inaccuracy  of  predicting  its  velocity.  Rather  than 
iterate  on  the  position  of  the  core,  the  analytical  predictions  of  Smith6 
were  used  for  the  positions  of  the  core  (Figure  22).  His  conical  flow  as- 
sumption is  followed  in  the  model  to  the  trailing  edge,  and  then  the  core 
is  assumed  to  trail  straight  back  with  a slight  incline,  as  in  the  fixed- 
wake  model . 

The  computational  procedure  for  the  free-wake  model  is: 

1.  Input  AR,  a,  Yy,  Zy,  Assumed  leading-edge  vortex 
strenqth  and  geometry.  Assumed  trai ling-edge  wake 
strength  and  geometry 
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2.  Calculate  and  decompose  the  wing-on-wing  influence  matrix 

3.  Calculate  vortex-  and  wake-on-wing  influence  matrices 

4.  Calculate  right-hand  side  of  Equation  23 

5.  Obtain  wing  doublet  distribution 

6.  Update  L.E.  vortex  and  T.E.  wake  strengths 

7.  Update  L.E.  vortex  and  T.E.  wake  positions 

(iterate  between  3 and  7 until  convergence) 

d.  Calculate  pressure  distribution  over  the  wing 
Some  exploratory  cases  run  as  described  indicated  that  a few  modifications  ' 

should  be  made.  The  first  modification  was  prompted  by  the  observation  that, 
when  one  of  the  free  vortices  was  of  odd  sign  in  strength,  geometry 
predictions  were  poor  to  the  point  of  being  unstable.  The  solution  was  to 
fit  the  strengths  of  the  strips  on  the  leading  :dge  sheet,  {K}y , to  a 
polynominal  curve.  By  forcing  these  strengths  to  fit  a curve  of  low  enough 
order  all  the  vortices  retained  the  same  sign.  In  fact,  for  a fourth  order 
curve,  the  vortices  were  all  of  the  same  sign;  and  yet  when  a converged  ca«e 
was  obtained,  the  leading-edge  boxes  and  attached  wake  strips  differed 
negligibly  ’n  strength,  implying  that  the  Kutta  condition  was  still  being 
satisfied.  A second  purpose  which  this  curve  fit  s '^ved  was  insuring  the 
doublet  strength  at  the  apex  was  zero. 

The  second  change  concerned  the  vortice..  as  they  leave  the  leading 
edge.  The  point  at  which  these  vortices  originate  is  on  the  leading 
edge  where  the  induced  velocities  cannot  be  accurately  calculated.  If 
the  tangent  off-flow  condition  is  arbitrarily  enforced,  that  Is,  the  net 
w velocity  at  the  leading  edge  set  to  zero,  the  lift  coefficient  will  change 
by  approximately  30  percent.  This  sensitivity  is  due  to  the  proximity  of 
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the  leading-edge  collocation  points  to  the  vortices  leaving  the  leading 
edge  and  has  been  noted  by  other  analytic  and  numerical  sources. 

• » 

The  modification  was  to  break  the  first  segment  of  each  vortex  into 
1 two  segments  with  the  break  occurring  at  20%  of  the  original  segment's 

length  (Figure  23).  By  breaking  the  segment  very  close  to  the  leading 
* edge,  another  point  is  established  at  which  the  velocities  can  be 

\ calculated  accurately,  which  leads  to  more  accurate  predictions  of  the  sheet 

geometry  in  this  critical  region. 

Convergence.  To  date,  no  objective  basis  has  been  established  for 
determining  when  the  program  has  converged  to  a solution.  The  possibility 
of  unrealistic  vortex  geometries  being  generated  still  dictates  that  a 
programmer  evaluate  the  program's  progress  every  five  to  ten  iterations, 

I as  opposed  to  approximately  twenty-five  iterations  required  for  convergence. 

Such  being  the  case,  it  has  been  left  to  the  progranmer  to  determine  convergence. 
This  determination  is  done  through  comparison  of  vortex  geometries  and  lift 
coefficients,  while  another  accurate  indicator  is  the  doublet  strength  at 
the  wing  tip.  Regardless,  judging  from  the  experience  of  having  run  the 
cases  presented  here,  further  iteration  should  vary  the  by  no  more  than 
5 percent. 

An  average  case  required  five  runs,  each  run  updating  the  geometry  of 
the  leading-edge  vortex  sheet  five  times.  The  total  run  time  of  75  seconds 
CPU  time  (on  the  Amdahl)  can  be  reduced,  s’nce  each  separate  run  requires 
decomposition  of  the  wing-on-wing  influence  matrix.  Another  factor 
affecting  computational  time  is  the  initial  approximation  to  the  sheet 
geometry.  Here,  results  were  first  obtained  at  a = 20°  for  aspect  ratios 


of  one  and  two,  using  hand-drawn  spirals  as  first  approximations. 
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Then  the  converged  geometries  of  these  cases  were  used  as  first  approximations 
at  a = 15°,  etc. 

Results.  The  results,  as  presented  in  Figures  24  throu  h 28,  are 
encouraging.  At  an  aspect  ratio  of  one,  the  vortex  lift  is  over  predicted 
(Figure  24),  but  the  general  trend  of  increasing  lift-slope  with  angle  of 
attack  is  followed.  At  aspect  ratio  two  (Figure  25),  the  comparison  with 
experiment  is  much  better.  The  difference  between  the  data  and  the  prediction 
is  within  the  uncertainty  band  of  the  data. 

With  a better  model  of  the  trailing  edge  the  center  of  pressure  (Figure  26) 
has  moved  forward  compared  to  the  fixed-wake  model  (Figure  15,  p.  30).  The 
free-wake  prediction  now  agrees  more  closely,  especially  in  trend,  with  ihe 
thick  wing  position,  but  still  does  not  show  the  thin  wing's  far  forward 
position. 

A carpet  plot  of  the  loading  coefficient.  Figure  27,  shows  the  predicted 
distribution  follows  the  peak  and  valley  form  of  experimental  measurement. 

Again,  as  in  the  fixed-wake  model,  the  loading  along  the  centerline  is 
predicted  well,  but  the  suction  peak  is  lower  and  further  outboard  of 
experiment.  This  discrepancy  between  the  positions  of  the  suction  peaks  may 
be  due  to  the  edge  shape  of  this  particular  wind-tunnel  model,  since  this  shape 
has  a noted  effect  on  the  core  position. 

Near  the  trailing  edge,  the  experimental  data  shows  a loss  of  vortex 
suction,  the  reason  for  the  forward  position  of  the  center  of  pressure  in 
thin  delta  wings.  This  loss  of  suction  seems  to  be  associated  with  the 
position  of  the  vortex  core.  As  It  nears  the  trailing  edge  the  core  passes 
out  of  the  conical  flow  region  and  starts  to  parallel  the  freestream. 
Consequently,  its  height  above  the  trailing  edge  is  greater  than  that 
predicted  by  conical  flow,  and  its  influence  is  therefore  less. 
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Smith's**  theory  fcr  predicting  the  core  position  should  not  be  over- 
looked as  a source  of  inaccuracy  in  the  pressure  distribution.  His  core 
positions  are  also  known  to  be  slightly  outboard  of  reported  data. 

Nevertheless,  some  aspects  of  the  theory  are  borne  out  by  experimental 
observation  and  niflnerical  check.  For  example,  conical  flow  predicts  the 
same  leading-edge  sheet  geometry  for  a unit  aspect  ratio  wing  at  a * 10° 
as  it  does  for  an  aspect  ratio  two  wing  at  20°.  Using  a fixed  sheet 
geometry  predicted  by  the  program  for  AR  - 2,  a = 20°,  the  unit  aspect 
ratio  pressure  distribution  was  computed.  The  lift  coefficient  obtained  this 
way  varied  by  less  than  5%  from  the  free-wake  case;  but,  the  error  in  the 
theoretical  position  of  the  core  seems  the  most  likely  cause  for  the  over 
prediction  of  the  vortex  lift  on  the  aspect  ratio  one  wing.  If  this  is  so, 
a simple,  fixed-wake  model  based  on  the  geometries  predicted  by  the  free-wake 
model  will  be  difficult  to  achieve  without  using  conical  flow. 

Figure  28  shows  one  example  of  the  predicted  geometry.  The  spiral  reflects 
more  of  the  experimentally  observed  roll-up  than  the  results  of  Kandil  et.  alJ® 
or  Rehbach.1^  A disadvantage  of  the  discrete  vortex  scheme  is  also  shown. 
Beneath  the  core  the  spanwise  velocities  are  so  high  that  no  more  than  one 
vortex  is  ever  present  between  the  core  and  the  wing.  Thus,  the  sheet  in  this 
area  is  not  well  represented;  which  may  be  the  cause  of  suction  peaks  lower 
than  observed. 

Future  Improvements.  The  free-wake  model  must  be  modified,  though, 
if  it  is  to  be  useful  for  higher  aspect  ratio  wings.  It  appears  that  for 
aspect  ratios  of  three  and  higher  the  calculations  for  vortex  segment 
; position  become  incorrect.  At  these  aspect  ratios  the  u velocity  beneath 

the  core  becomes  negative.  The  equations  (25,  26)  for  predicting  the 
| position  of  the  segments  implicitly  assume  positive  u velocities  and,  hence. 
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they  tail  at  higher  aspect  ratios.  A new  coordinate  system  using  the 
vortex  core  as  one  axis  and  radial  and  angular  coordinates  similar  to  those 
described  will  probably  be  adopted.  The  velocity  along  the  core  is  known 
to  be  greater  than  the  freestream  which  will  alleviate  the  problem  of  negative 
u velocities. 

At  the  same  time  it  would  be  advantageous  to  modify  the  grid  on  the 
wing,  since  at  low  values  of  a,  (a/ tan  e),  the  wing  doublet  strength  beneath 
the  core  has  a very  high  gradient.  Such  gradients  are  contradictory  to 
the  assumption  that  K may  be  considered  constant  in  each  panel  on  the  wing. 

Therefore,  a fine  grid  should  be  imposed  on  the  wing  beneatr.  the  core,  which 
would  also  assist  in  calculating  the  loading  in  this  critical  area. 

A more  sophisticated  predictor  will  allow  the  program  to  reach  convergence 
without  operator  help.  This  will  improve  run  mes;  but  it  will  require 
the  adoption  of  convergence  criteria.  The  use  of  convergence  criteria  will 
also  be  advantageous  since  the  introduction  of  objectivity  could  improve 
the  consistency  of  the  program. 

Finally,  the  vortex  core  should  be  positioned  on  a basis  other  than 
conical  flow  theory.  Experimental  results  might  be  used,  but  ultimately 
a numerical  prediction  will  be  made. 

Conclusions 

A general  potential  flow  theory,  including  the  effects  of  leading- 
edge  separation,  has  been  presented  and  applied  specifically  to  delta  wings. 

Two  different  models  of  the  leading-edge  separation  were  constructed  numeri- 
cally: a fixed-wake  model,  and  a more  refined  free-wake  model. 

The  fixed-wake  model  is  simple,  fast,  and  accurate  for  low  aspect  ratio 
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wings,  but  it  does  not  provide  detailed  information  about  the  flow. 
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The  free-wake  model  has  a larger  aspect  ratio  range  for  which  compara- 
tively accurate  answers  can  be  obtained.  It  also  provides  some  qualitative 
results  on  the  shape  of  the  leading-edge  vortex  sheet.  Its  disadvantages 
are  longer  computational  times  and  the  need  for  operator  supervision. 

Methods  of  enlarging  the  range  of  aspect  ratios  for  which  the  models 
are  valid  have  been  proposed.  Some  of  these  are  now  being  evaluated  and 
will  be  incorporated  into  a generalized  program  combining  the  fixed-wake 
and  free-wake  concepts. 
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II.  EXPERIMENTAL  INVESTIGATION 
Equipment  and  Tests 


Most  of  the  work  on  delta  wings  in  unsteady  flow  has  been  done  by  moving 
(33-421 

wings  in  a steady  flow'  '.  Whilst  the  basic  mechanics  of  this  approach 
are  straightforward  it  does  lead  to  problems  with  data  aquisition  and  recording 
and  might  not  be  representative  of  an  aircraft  entering  a gust.  In  particular, 
the  thin  wing  used  in  these  tests  was  subject  to  unacceptable  flexure  when 
pitched  rapidly.  Therefore,  the  reverse  approach  of  statically  mounting 
the  wing  in  a variable  airstream  was  pursued. 

For  the  current  tests,  a thin  sharp-edged  delta  wing,  AR  = 2 , was 
mounted  in  the  Texas  A&M  University  7'  x 10'  Low  Speed  Wind  Tunnel  (Fig.  29) 
and  a stream  oscillation  devi ce^ 43-44 ^ (Fig.  30)  used  to  produce  a sinusoidal 
variation  in  flow  angle  of  the  airstream.  The  system  uses  two  flapped  vanes 
at  the  entrance  to  the  test  section,  the  sides  of  which  have  been  removed. 

The  flaps  are  oscillated  by  a variable  speed  electric  motor  and  produce  a 
sinusoidally  oscillating  airstream,  of  constant  amplitude,  in  the  test  section. 

Frequencies  of  oscillation  of  2Hz  and  4Hz  were  used  and  in  conjunction 
with  an  airspeed  of  83.8  ft. /sec  resulted  in  reduced  frequency  parameters, 
k,  of  0.2  and  0.4  at  a Reynolds  number  of  1.3  x 10^.  With  these  frequencies 
and  tunnel  speed  the  total  wavelengths  of  the  flow  were  42  ft  and  21  ft  as 
compared  to  the  root  chord  length  of  4 ft  for  the  delta  wing.  The  flow 
angle  amplitude  chosen  for  the  tests  was  +_  8°  about  a mean  angle  of  attack 
of  15°. 

The  tests  were  divided  into  two  sections,  flow  visualization  and  pressure 

measurements.  For  the  flow  visualization  a liauid  nitrogen  and  steam  smoke 

(45) 

generator  developed  at  Texas  A&M'  ’ was  used  to  seed  the  vortex  cores. 
Initially  the  position  of  the  vortex  core  and  burst  point  (if  it  occurred) 
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in  steady  flow  were  determined  at  5°,  10°,  15°,  20°,  25°,  30°  angle  of 
attack,  by  photography  (Fig.  31).  In  the  unsteady  case  photographs  of 
the  vortex  locations  at  a particular  point  in  a cycle  were  obtained  using 
strobe  lights,  slaved  off  the  flap  drive  motor,  to  illuniwate  the  flow 
(Figs.  32-34)  the  rest  of  the  tunnel  being  darkened.  The  loci  of  the  vortex 
burst  point  versus  angle  of  attack  were  determined  from  these  photographs 
(Fig.  35).  For  detailed  chronological  analyses,  the  two  unsteady  tests 
were  filmed  using  high-speed  cameras  opp “icing  at  200  frames  per  second. 
These  f’lms  can  be  correlated  with  other  data  by  a flap-motor  indicator 
(reflected  by  a mirror)  in  the  camera's  view. 

For  the  pressure  measurements  twenty  Validyne  DP9  pressure  transducers 
were  mounted  on  the  wing  with  one  side  of  each  transducer  connected  to  a 
port  on  the  upper  surface  (Fig.  36)  while  the  other  side  was  connected  to 
a plenum  referenced  to  free-streair  static  pressure.  Outputs  from  all 
transducers,  together  with  a signal  indicating  the  flap  position  were 
recorded  on  Honeywell  visicorders.  Again  results  were  obtained  for  both 
steady  and  unsteady  flow.  After  reduction  the  load  distributions  on  the 
wing  were  plotted  at  several  angles  of  attack  throughout  a typical  cycle 
at  each  frequency  (Fig.  37). 


Discussion 

1)  Flow  Visualization 
a)  Steady  Flow 

At  5°  angle  of  attack  in  steady  flow  (Fig.  31)  the  vortices  above 
the  wing  are  weak  and  it  proved  impossible  to  seed  the  cores  adequately 
with  smoke.  When  a is  increased  to  10°  the  vortex  core  is  well  defined 
and  easily  seeded,  towards  the  trailing  edge  the  core  tends  to  curve 


into  the  stream  direction.  Increasing  , to  16°  causes  the  burst  to  move 
jr  to  the  trailing  edge  in  the  lower  vortex,  no  bursting  being  evident  ir 
the  'jo per  vortex.  The  slight  i ymne*r>  is  probably  due  to  either  a small 
yaw  angle  on  the  wing  or  a disturbance  to  the  *low  caused  by  the  smoke  probe 

At  20°  angle  of  attack,  the  vortices  have  moved  further  inboard  ana 
the  burst  point  nas  moves  *orward  to  almost  SO3,  of  tne  root  chord.  These 
trends  continue  until  at  j = 30°  the  burst  point  has  moved  fo-ward  to 
almost  tne  2Q"  root  chord  position. 

Tne  position  of  tne  burst  point  as  a function  ,«f  angle  of  auack  i: 
plotted  in  figure  35  ana  as  can  be  seen  the  burst  point  moves  stadily 
•cor\:ard  with  increasing  angle.  The  points  at  which  the  :ore  begins  to  burst 
and  when  the  smoke  is  no  longer  visible  are  plotted  to  indicate  the  severity 
of  the  burst.  Numerically,  the  results  are  cons  .tent  with  Lambourne  s . at 
all  angles  of  attack  the  burst  point  in  these  tests  lies  slightly  anead  cf 
Lamboume ' s^J  measurements  as  one  would  expect  for  a slightly  larger  aspect 
ratio  (2.0  as  opposed  to  1.86). 

In  the  current  steady  flow  tests,  the  burst  was  always  of  the  axi sym- 
metric "tulip"  type  (this  was  also  true  for  some  preliminary,  tests  at 
5 IF 

- 3 x 10  ).  According  to  Sarpkaya  ' tie  spiral  type  of  vortex  breakdown 
is  more  commonly  observed  over  delta  wings  _id  Maltbv  et.al.  £ were  able  to 
produce  a spiral  burst  over  an  aspect  ratio  one  wing  in  a wind  tunnel  at  Rey- 
nolds numbers  of  1.5  and  6 x 10^.  However  the  observations  reported  nere  art- 
consistent  with  tne  results  of  Lambourne  who  pointed  out  mat  his  observations 
of  the  spiral  type  of  burst  were  at  very  low  Reynolds  numbers  in  a wate*-  tume 

b)  Oscillatory  Flow 

The  analysis  of  the  unsteady  flow  is  based  t the  still  photos  presented 
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in  Figures  32  and33and  the  high  speed  films.  The  angle  of  attack  is 
referred  to  the  instantaneous  angle  of  attack  at  the  aerodynamic  renter 

( »2/3c  = + s’n 

The  still  photos  of  the  flow  at  2Hz  (Fig.  32)  indicate  that  at  15.3° 
with  the  angle  of  attack  increasing  the  vortex  is  still  fairly  weak  and 
no  burst  occurs  above  the  wing.  When  the  angle  has  increased  to  20.9° 
the  vortex  has  moved  inboard  and  bursts  at  about  the  60t  root  chord  point. 

It  should  be  noted  however  that  this  burst  does  not  appear  to  be  of  the  same 
type  as  that  observed  in  steady  flow  (Fig.  31),  the  burst  does  not  have  the 
classic  'tulip'  shaped  bubble  normally  seen  with  an  axisymmetric  burst. 

Rather  the  core  seems  to  expand  gradually,  become  turbulent  and  dissipate. 

By  the  time  the  maximum  angle  of  attack  has  been  reached  23°  the  burst 
point  has  moved  forward  but  still  appears  to  be  of  a gradual  diffusion  type. 

At  20.5°  and  decreasing  the  burst  point  has  not  changed  location  but  its 
character  has  changed.  This  trend  continues  until  a = 14.7°  when  the  burst 
has  moved  downstream  and  become  a helical  type  burst  (see  also  Fig.  34)  where 
the  core  has  an  abrupt  kink  and  then  swirls  around  its  axis  before  dissipating 
into  turbulence.  The  whole  process  is  much  more  rapid  than  the  gradual 
increase  of  core  size  previously  described.  By  9.1°  angle  of  attack  the 
burst  point  is  near  the  trailing  edge  of  the  wing  and  the  vortex  is  decreasing 
in  strength.  As  the  angle  passes  through  a minimum  the  burst  moves  off  the 
wing  and  the  core  moves  outboard  before  moving  inboard  again  with  increasing 
angle  (i  - 9.5°). 

The  point  of  the  burst  for  k = 0.2  significantly  lags  the  position  in 
steady  flow,  also  the  burst  point  in  oscillatory  flow  never  moves  as  forward 
on  the  wing  as  it  does  in  steady  flow.  Tnis  is  presumably  because  it  takes 
a finite  time  for  the  vortex  strength  to  increase  after  the  leading  edge 
shedding  rate  has  changed  due  to  angle  of  attack  increase,  it  will  then  take 


a further  period  of  time  for  the  instabilities  to  cause  the  burst  point 
to  move  forward.  By  this  time  the  angle  of  attack  has  decreased  and  the 
shedding  rate  changed,  i.e.  the  wing  is  not  at  the  maximum  angle  of  attack 
long  enough  for  the  steady  state  conditions  to  occur,  the  integrated  vortex 
strength  which  must  affect  the  burst  location  is  always  less  in  unsteady 
flow  than  it  is  at  the  maximum  angle  of  attack  in  steady  flow. 

Figure  38  produced  from  data  taken  from  the  high  speed  film  presents 
graphically  the  position  and  character  of  the  burst.  Firstly,  it  should 
be  noted  that  the  vortex  core  was  not  seeded  well  enough  to  allow  plotting 
of  the  burst  during  the  increasing  part  of  the  cycle.  Up  to  the  maximum 
angle  the  diffuse  core  focuses  until  the  beginning  and  end  of  the  burst 
can  be  determined  with  some  certainty.  At  #'^°  the  core  begins  bursting 
at  601  of  the  root  chord.  The  initial  burst  point  moves  only  slightly 
forward  while  the  terminal  end  of  the  burst  moves  rapidly  forward  as  the 
burst  changes  from  gradual  diffusion  to  a helical  burst.  As  the  angle  of 
attack  decreases  from  17.5°  the  burst  moves  rearward  until  the  core  is  too 
weak  to  appear. 

For  k = 0.4  f^ig.  33)  the  trends  are  similar  to  those  observed  at  k = 0. 
At  17.3°  with  the  angle  of  attack  increasing  there  is  no  burst  above  the 
wing.  As  the  angle  increases  to  just  above  20°  the  burst  moves  onto  the 
wing  and  by  22. 1°  it  occurs  at  about  the  601  chord  location.  Determination 
of  the  exact  burst  point  at  this  angle  is  extremely  difficult  because  of 
the  very  gradual  expansion  of  the  core  into  turbulence.  Just  after  the 
maximum  angle  has  been  reached  (a  = 22.7)  the  burst  point  is  much  better 
defined  but  is  still  axi symmetric  in  nature.  As  the  angle  decreases 
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through  18.8°,  12.7°  and  7.9°  the  burst  changes  character  becoming  helical 
and  moving  downstream.  The  burst  point  moves  off  the  wing  as  the  minimum 
angle  of  attack  is  reached. 

Figure  39  is  similar  to  Figure  38  but  the  diffuse  phase  of  the  vortex 
is  much  more  exaggerated.  Again,  the  determination  of  the  initial  and 
terminal  burst  points  is  difficult  near  the  maximum  angle.  Presumably 
the  effective  angle  at  the  apex  increases  quickly  enough  that  the  initial 
burst  does  not  propagate  upstream  but  is  reformed  each  cycle  and  moves 
rearward  with  decreasing  angle  of  attack. 

2)  Pressure  Measurements 

All  the  pressure  data  obtained,  both  steady  and  unsteady,  is  presented 
in  the  form  of  carpet  plots  at  several  angles  of  attack  in  Fig.  37. 

With  the  wing  at  15°  angle  of  attack  in  steady  flow  (Fig.  37a)  the 
two  suction  peaks,  near  the  nose  are  narrow  and  close  to  the  leading  edges, 
further  downstream  the  peaks  move  slightly  inboard,  spread  out  and  decrease 
in  magnitude  due  to  the  effects  of  secondary  separation  and  the  Kutta 
condition  of  the  trailing  edge.  In  unsteady  flow  the  distribution  is 
similar  but  the  peaks  are  slightly  lower  for  k = 0.2  and  lower  still  for 
the  k = 0.4  case.  This  trend  continues  up  through  a = 17.5°(Fig.  37b  - Note 
steady  flow  data  was  only  obtained  at  5°  intervals)  with  the  suction  peaks 
being  higher  at  k = 0.2  than  at  k = 0.4. 

At  a = 20°  (Fig.  37c)  these  trends  are  still  true  at  the  most  forward 
measuring  station  (x/c^  = 0.4)  but  at  the  midstation  (x/c^  = 0.6)  the  peaks 
in  unsteady  flow  are  slightly  greater  than  in  steady  flow.  The  -eason  for 
this  is  probably  in  the  vortex  burst,  in  steady  flow  the  burst  occurs 
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ahead  of  whilst  In  unsteady  flow  it  has  not  yet  reached,  this  location, 
thus  in  unsteady  flow  full  vortex  suction  is  still  being  maintained. 

Peaks  for  k ■ 0.2  are  slightly  higher  than  for  k = 0.4.  At  the  rear 
station  the  peaks  for  k = 0.2  and  0.4  are  about  the  same  magnitude  as 
the  vortex  has  burst  upstream  of  this  location  at  k = 0.2  but  has  not  for 
k = 0.4,  thus  the  loss  of  suction  due  to  bursting  cancels  the  loss  due 
to  phase  lag. 

By  the  time  maximum  angle  of  attack  (23°Fig.  37d)  has  be^n  reached 
the  vortex  burst  has  moved  ahead  of  x/c^  =0.6  for  bo. h frequency  parameters 
and  distributions  of  x/cr  = 0.6  and  0.8  look  similar  for  both  cases. 

The  suction  peaks  at  x/c^  = 0.4,  where  no  burst  occurs  in  either  case,  are 
greater  for  the  lower  frequency. 

By  comparing  the  distributions  at  the  various  angles  of  attack  it  can 
be  seen  that  the  total  load  on  the  wing  increases  steadily  as  the  angle 
of  attack  increases  from  15°  to  23°. 

When  a decreases  to  20°  (Fig.  37e)  the  loading  for  all  cases  is 
similar,  and  might  be  expected  since  this  is  the  angle  at  which  the  burst 
patterns  are  in  phase.  A major  difference  between  the  steady  and  unsteady 
data  is  that  the  dips  due  to  secondary  separation  apparent  in  the  steady 
data  are  not  present  in  the  unsteady  distributions,  this  might  just  be  due 
to  the  ports  being  located  at  the  'wrong'  places  to  measure  secondary 
separation  in  unsteady  flow.  At  the  most  forward  station  the  peaks  for  k = 0.2 
lie  slightly  outboard  and  are  slightly  lower  than  the  steady  or  k = 0.4 
results.  The  distribution  changed  little  between  23°  and  20°. 

Again  little  change  in  loading  occurs  as  the  angle  of  attack  further 
decreases  to  17.5°  (Fig.  37f)  but  as  the  angle  decreases  to  15°  (Fig.  37g) 


the  load  on  the  wing  begins  to  decrease.  An  interesting  point  at  this 
angle  is  the  distributions  of  x/cR  * 0.6  and  0.8,  the  suction  peaks  In 
unsteady  flow  are  quite  broad  whilst  they  are  sharp  and  narrow  for  steady 
; flow.  In  steady  flow  at  this  angle  vortex  bursting  does  not  occur  until 

the  trailing  edge  whilst  at  k = 0.2  it  occurs  at  about  x/cR  * 0.6  and 
somewhat  further  forward  for  k = 0.4  causing  the  peaks  to  become  broader 
and  lower.  A similar  phenomena  occurs  at  a * 12.5°  (Fig.  37h)  where  the 
; peaks  at  x/cR  =0.6  and  0.8  are  broader  and  lower  for  the  higher  frequency. 

| As  a decreases  to  iO°(Fig.  37i)  the  total  lift  continues  to  decrease. 

I 

i The  vortex  loading  in  unsteady  flow  at  the  two  forward  stations  is  considerably 

| greater  than  for  steady  flow.  At  the  rear  station  vortex  burst  causes 

I the  peaks  to  broaden  and  decrease  in  magnitude  for  unsteady  flow. 

t At  the  minimum  angle  of  attack,  a = 7°  (Fig.  37j),  the  vortex  peaks 

at  x/cR  = 0.4  for  k = 0.2  have  become  very  narrow  and  have  not  moved  as 

far  outboard  as  have  the  peaks  for  k = 0.4.  The  reason  for  this  is  not 

known  but  it  mignt  be  connected  with  the  rapid  collapse  of  the  vortex  core 
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at  low  angles  in  large  amplitude  oscillatory  flow  as  discussed  by  Lowson 
As  the  angle  of  attack  increases  back  up  to  10°  (Fig.  37|0  the  lift 
continues  to  decrease.  It  is  not  until  12.5°  (Fig.  37t)  that  the  lift 
starts  to  increase  again  with  increasing  magnitudes  of  vortex  suction. 

In  general  it  can  be  seen  that  in  oscillatory  flow  there  is  a phase 
lag,  that  increases  with  frequency,  in  the  development  of  vortex  lift. 

This  is  not  so  in  the  loading  away  from  the  influence  of  the  vortices  i.e. 
along  the  center  of  the  wing.  Flow  that  comes  over  the  vortices  and  reattaches 
along  the  centerline  of  the  wing  (the  'potential'  flow)  responds  fast  enough 
that  there  are  no  phase  lags  within  the  accuracy  of  measurement  'except  at 

» 
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high  angles  where  the  center  section  is  influenced  by  the  vortices).  The 
phase  lag  that  occurs  in  the  vortex  loading  must  occur  because  of  the 
finite  time  it  takes  the  vortex  strength  to  change  after  a change  in  the 
leading  edge  shedding  rate,  this  is  consistent  with  the  measured  lags  in 
the  locus  of  the  burst  points. 

In  the  results  presented  here  no  wind  tunnel  boundary  connections  have 
been  applied.  Estimates  for  the  magnitude  of  the  angle  of  attack  correction 
were  made  using  standard  steady  flow  vortices  of  lift  coefficient  and  these 
indicated  that  at  a = 7°  and  23°  correction  of  erder  -0.2°  and  -0.6°  would 
De  needed.  However  this  type  of  corrections  cannot  be  rigorously  justified 
in  unsteady  flow  and  so  none  have  been  applied. 

Conclusions 

In  steady  flow  at  Reynolds  numbers  of  order  1 x 106,  a delta  wing  of 
AR  = 2 the  vortex  bursting  that  occurs  is  of  the  axisymmetric  'tulip'  type 
and  the  burst  point  moves  upstream  with  increasing  angle  of  attack.  In 
unsteady  fow  = 15°  + 8 sin  wt)  this  type  of  burst  does  not  occur.  As  the 
angle  of  attack  increases  the  vortex  core  appears  to  undergo  a gradual 
viscous  dissipation  into  turbulence  whilst  with  decreasing  angle  a helical 
(i.e.  the  vortex  core  swirls  about  its  own  axis)  burst  was  observed.  The 
locus  of  the  burst  point  in  unsteady  flow  lags  the  steady  position  over  most 
of  the  cycle.  The  burst  point  Joes  not  move  as  far  forward  in  unsteady  as 
it  does  in  steady  flow,  probably  because  the  vortex  strength  takes  a finite 
time  to  build  up  and  the  leading  edge  vorticity  shedding  rate  is  not 
maintained  at  the  maximum  value  long  enough  for  tne  vortex  to  reach  equilibrium. 


Measured  pressure  distributions  indicate  that  the  suction  pressures 
due  to  the  vortices  in  unsteady  flow  lag  behind  the  steady  state  pressures 
through  most  of  the  cycle  whilst  the  loading  in  regions  of  'potential' 
flow  remain  in  phase  with  the  angle  of  attack  variation  confirming  that 
the  vortices  take  a finite  time  to  change  strength  after  a change  in  angle 
of  attack.  Vortex  bursting  decreases  the  suction  pressure  and  broadens  the 
peak  causing  a decrease  in  lift  downstream  of  the  burst.  When  the  angle  of 
attack  is  decreasing  the  locus  of  vortex  bursting  lags  the  steady  state 
position  by  large  amounts  causing  a significant  loss  of  lift  over  the  rear 
of  the  wing  to  much  lower  angles  of  attack  then  might  be  anticipated  from 
steady  flow  data. 
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